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Abstract 

We discuss the imphcations for cosmic microwave background (CMB) observables, of a class of 
pre-inflationary dynamics suggested by string models where SUSY is broken due to the presence 
of D-branes and orientifolds preserving incompatible portions of it. In these models the would-be 
inflaton is forced to emerge from the initial singularity climbing up a mild exponential potential, 
until it bounces against a steep exponential potential of "brane SUSY breaking" scenarios, and 
as a result the ensuing descent gives rise to an inflationary epoch that begins when the system is 
still well off its eventual attractor. If a pre-inflationary climbing phase of this type had occurred 
within 6-7 e-folds of the horizon exit for the largest observable wavelengths, displacement off 
the attractor and initial-state effects would conspire to suppress power in the primordial scalar 
spectrum, enhancing it in the tensor spectrum and typically superposing oscillations on both. 
We investigate these imprints on CMB observables over a range of parameters, examine their 
statistical significance, and provide a semi-analytic rationale for our results. It is tempting to 
ascribe at least part of the large-angle anomalies in the CMB to pre-inflationary dynamics of 
this type. 
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1 Introduction 



An epoch of primordial inflation [Ij is widely accepted to have premised the Big Bang, since 
it readily accounts for the flatness, homogeneity, isotropy and causal structure of the initially 
thermalized plasma, as well as for the apparent absence of dangerous topological relics in our 
Universe. In its simplest realizations, inflation is also remarkably successful in accounting for 
the observed CMB anisotropies, and provides a quantum mechanical origin for the primordial 
inhomogeneities around which we believe that all visible structure formed. 

In order to solve the flatness problem, inflation must have lasted for a certain number of 
e-folds, with a lower bound that depends on the putative reheating temperature of the Universe 
(at least 58 e-folds for GUT-scale reheating). In toy models where the inflaton is treated as 
a classical scalar field minimally coupled to gravity, substantially more than this minimal e- 
folding can readily be attained with fine-tuned potentials or large field excursions [2]. Matters 
are markedly different, however, if one demands that infiation result from an effective field theory 
(EFT) where the infiaton is an effectively light degree of freedom, since loop corrections tend 
to make long periods of inflation unnatural. If the effective potential, rather than the classical 
potential, is to be responsible for the phenomenon [3], various tunings are in fact to be made 
or various symmetries are to be invoked, which are furthermore required to survive the UV 
completion of the theory. One of the most challenging issues for model builders is perhaps the 
"eta problem", wherein corrections to the inflaton mass of order H^, where H is the Hubble 
parameter, are readily generated by operators of dimension six in the EFT, while the inflaton 
mass scale ought to be far lower for successful inflation (see |4j for a discussion of the problem 
and a proposal to evade it). In spite of many suggestions to bypass it in specific contexts (see 
also [5]), the eta problem seems to persist in generality [6^. 

Little is known about the actual EFT in which infiation is to be embedded, and as we have 
stressed, sustaining it for the needed duration appears to be a delicate proposition in realistic 
constructions. In low-energy supergravity models inspired by String Theory, the challenge is 
to arrange for the infiaton to undergo slow roll in a region of the potential that can sustain 
it without overshooting [8j , while in general EFT treatments [3l [9] the challenge is to arrange 
for effective potentials leading to infiaton masses that are much smaller than the Hubble scale. 
All these requirements are ever more challenging if one requires direct couplings to a matter 
sector that contains the Standard Model, which is needed if reheating is to occur efficientl}{l. In 
light of these considerations, it behooves us (as it has behooved other authors [TT]) to entertain 
the possibility that infiation did not last that much longer than is needed to solve the fiatness 
problem. 

Long-wavelength distortions of CMB observables are obvious corollaries of short infiation. 
Although we are limited by cosmic variance at the largest angular scales, some statistically signif- 
icant imprints of a possible pre-inflationary phase could in fact persist in the sky if infiation did 
not last more than 6-7 e-folds longer than it had to. For instance, a single pre-infiationary par- 
ticle could have sourced ring-like structures in the CMB [12^ II 3 j^ . Alternatively, a non- vanishing 
CMB dipole sourced by super-horizon inhomogeneities that were not completely washed away, 
although insignificant for the CMB power spectrum, could have had sizable consequences for 

^See also [7] for a discussion of the eta problem as it relates to uplifting in supergravity theories. 
See |10j for interesting observational consequences that could result in the context of large volume compacti- 
fications. 

''The many-particle superposition of this effect would bring about an excess of power (i.e. a feature) at long 
wavelengths in the power spectrum. 
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the large-scale structure of the Universe. To wit, there may be evidence for a large-scale bulk 
flow [H] that is in significant discrepancy with the predictions of ACDM models, which could 
readily be explained if the CMB possessed an intrinsic dipole moment |15] . Taken together 
with statistically significant observations |16j of the lack of power at large angles in the angular 
correlation function (as opposed to the angular power spectrum j^, this suggests that we may be 
seeing something other than the vanilla initial conditions generated by a phase of long inflation. 

With this in mind we explore here the consequences for CMB observables, of a genre of 
pre-inflationary dynamics |17j that presents itself in String Theory [18], in orientifold vacua 
|19j where supersymmetry (SUSY) is broken by classically stable combinations of (anti) D- 
branes and orientifolds [2Ql- The peculiarity of these models is that a scalar field is typically 
forced to emerge from the initial singularity climbing up an exponential potential, so that it can 
subsequently inject a slow-roll phase that begins off the eventual attr actor [17] [j. As we shall 
see, aside from providing a natural mechanism to initially displace the infiaton up its potential 
(so that the origin of its initial potential energy is accounted for), this climbing phase and the 
transients associated with it conspire to modify the power spectrum at long wavelengths. We 
find suppression of power at the largest scales for the scalar spectrum, (mild) enhancement 
for the tensor spectrum, with superposed oscillations at long wavelengths in both cases. We 
anticipate that a better fit to CMB observations may result and account at least for part of the 
large-angle anomaly previously discussed. 

The relevant mechanism, originally termed "brane SUSY breaking" [20], results from the 
presence of (anti)branes and orientifolds preserving incompatible fractions of the original SUSY, 
which is dictated by charge constraints in some orientifold vacua of String Theory. This gives 
rise to tree-level exponential potentials, whose lack of critical points would bring about in 
a flat background, tadpoles signalling the need for complicated resummations [271 EH]- In a 
cosmological setting, however, the same potentials have the intriguing consequence that a scalar 
fleld can be forced to emerge from the initial singularity while climbing up. As discussed in |17j . 
this is already the case in the simplest instance of "brane SUSY breaking" , the ten-dimensional 
Sugimoto model [29j, where the initial climbing can even inject an eventual slow-roll phase 
driven by a milder exponential [3^ of non-perturbative origin. Similar phenomena take place 
in the richer four-dimensional KKLT scenario [31], whose "uplift" may also be ascribed to the 
phenomenon of "brane SUSY breaking" and gives rise again to the climbing behavior. From this 
viewpoint, KKLT-like scenarios are also interesting since, as shown in [ITJ, the combined effects 
of a climbing phase and cosmological friction could have trapped some moduli in metastable 
vacua right after the initial singularity. A natural objection to this type of analysis is that 
even restricting the attention to tree-level potentials, which are expected to dominate at weak 
coupling in String Theory, curvature corrections near the initial singularity could invalidate the 
climbing mechanism. We are well aware of the problem, and we do not have as yet a truly 
satisfactory answer to this important question. Still, we are persuaded that a mechanism that is 
capable of naturally injecting inflation deserves to be investigated further, and it is conceivable 
that our setup will survive the relevant corrections at least in special classes of models. Moreover, 
the detailed consequences for CMB observables that we are about to describe depend largely 
upon the dynamics of the climbing inflaton at later times, when it idles about its eventual 

''Taking into account the observed alignments at low multipoles, the chances that such large-angle anomalies 
would be produced by a concordance model has been estimated by the authors of [16] as 1 in 10^. 

^We thank E. Kiritsis for calling to our attention some previous work on exponential potentials when the 
results of [17] were presented at "String Phenomenology 2010". A vast previous literature actually includes 
[211 1221 1231 1241 1251 126) and is mostly devoted to the asymptotic behavior of similar classes of models even with 
more scalar fields, but Russo discussed in [26] the classical solutions that we about to review in Section [J] 
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attractor. Although clearly embodied by specific string-inspired constructions where a single 
minimally-coupled scalar field is subject to the relatively simple class of potentials 

V{^) = ^ (e-^<^ + e-^^^) , (1.1) 

where a "hard" exponential brought about by brane SUSY breaking accompanies a slow-roll 
term (with < 7 < l/\/3) that dominates the eventual inflationary dynamics, this type of 
phenomenon ought to be relevant in a wider context. As we shall see, the scalar is forced to 
emerge from the initial singularity climbing up the mild portion of the potential, so that at some 
point it experiences the hard exponential of String Theory, essentially bouncing against it. The 
ensuing dynamics equilibrates rather slowly to the eventual attractor, and this brings about a 
sizable infrared suppression of the CMB power spectrum. Our result is not completely new, 
since an infrared suppression was also related in [32] to a pre-inflationary fast roll in double-well 
potentials. As we shall see, however, while a suppression is generic in models with an initial 
singularity, it is particularly sizable for these scalars that are initially forced to climb until they 
bounce off a hard exponential "wall", with potentially observable consequences. 

We can summarize our findings as having uncovered: 

• a statistically significant suppression of angular power over a range of lower multipoles, 

• which reflects pre-inflationary dynamics that settled onto slow-roll up to several e-folds 
before the present Hubble scale exited the horizon. 

Our construction was motivated by a class of string constructions: 

• where the breaking of SUSY is closely tied to the string scale; 

• which provide naturally the very impetus for inflation to begin; 

• which motivate a new class of phenomena - starting inflation off the attractor. 

The paper opens in Section [2] with a brief review of how a single minimally-coupled scalar field 
with an exponential potential of large enough slope is forced to emerge from the initial singularity 
while climbing up. We leave out a discussion of the more complicated KKLT-type constructions 
(for which the reader is referred to ^7\), since the ensuing analysis of CMB perturbations is 
restricted, for the sake of simplicity and predictivity, to models whose cosmological evolution 
is dominated by a single scalar field. A closer look at the dynamics of this scalar field, which 
is presented in the following sections paying attention to both pre-inflationary and slow-roll 
regimes, shows that inflation can effectively begin long before the solution comes close to the 
eventual attractor. In Section [3] we begin to analyze the behavior of correlation functions in 
this class of backgrounds, starting from a perturbative treatment of the corrections to CMB 
observables induced by inflating off the attractor, and in Section [J] we review some key properties 
of scalar and tensor power spectra. In Section [5] we numerically evaluate these power spectra for 
various choices of parameters of the climbing/inflationary potential, present an understanding 
of our results in terms of spectral running, and discuss the prospects for offering a better flt to 
observations. In Section [6] we elaborate upon a general, albeit semi-quantitative, WKB picture 
of the observable consequences of a "hard" exponential and the associated climbing phase, 
relying on some similarities between the Mukhanov-Sasaki and Schrodinger problems. Section [7] 
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contains our conclusions. Finally, Appendix [A] summarizes general features of cosmological 
perturbations, Appendix [B] describes an analytic model for scalar and tensor power spectra and 
Appendix [C] addresses the issue of initial-state effects on their long-wavelength tails. In the 
paper we shall refer to three types of time derivatives, with respect to a parametric time r that 
makes it possible to obtain some useful analytic information on the system, with respect to the 
cosmological time t and with respect to the conformal time rj. Following a standard practice 
in Cosmology, the latter two will denoted by "dots" and "primes" respectively, while the first 
type of derivatives will be indicated explicitly. These conventions differ from those of [l7j, where 
"dots" referred to r-derivatives. 



2 Climbing scalars and String Theory 



Let us begin with a brief review of some results of [17] • Our analysis rests on a class of low-energy 
four-dimensional actions of the typ^ 



S 



d X yj—g 



(2.1) 



for a single scalar field coupled to gravity, where for notational simplicity we readily turn to 
Planck units, setting from now on = SttGn = 1. One can study an interesting class of 
cosmological solutions letting 



and making the convenient gauge choice [22^ ^^ [26] 



where M is a mass scale related to the potential V{cl)) 0. With the redefinitions 



— \l n 'P 1 A = "iA , 
and confining our attention to an expanding Universe, so that 



one is thus led to recast the usual dynamical equation 

'(f) + ?,H^ + 1/^ = 0, 
where "dots" denote derivatives with respect to the cosmological time and 



H 



-L,/2y + - 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



In [T7| we were abiding to standard conventions in String Theory. Here we change slightly the Lagrangian, 
in order to work with a canonically normalized scalar field. 

^In general units, this gauge choice would involve the dimension-four ratio 
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into the convenient form 




) 



2 



+ 




) 



2 




(2.8) 



that only involves the logarithmic derivative of the potential V, and thus anticipates the key 
role of exponential potentials for these systems. 

Potentials of the type 



depending on a pair of exponential functions, with < 7 < 1 and where a positive relative 
coefficient could be absorbed by a shift of if, are naturally present in String Theory [18j and 
particularly so in orientifold models [19] with "brane SUSY breaking" |20j. In these models 
SUSY, broken at the string scale because charge conservation requires that (anti) D-branes and 
orientifold planes respecting incompatible portions of it be present, is non-linear ly realized a la 
Volkov-Akulov [33J in the low-energy supergravity [33]. In flat space exponential potentials of 
this type are somehow a nuisance, since their lack of critical points brings about tadpoles and 
thus the need for complicated vacuum redefinitions [271 [28] , but in cosmological settings they 
give rise to some interesting effects that we are about to recall. 

This type of scenario finds its simplest realization in the ten-dimensional Sugimoto model [29], 
where ip can be related to the dilaton. In this case the first term in ()2.9p reflects (projective) disk 
contributions accounting for the D9 anti-branes and the 09 planes present in the perturbative 
vacuum while the second accounts for the non-BPS D3 brane that was identified in [30] following 
|35j . The setting is also relevant, however, for the four-dimensional KKLT scenario [31], where 
as stressed in [17] F-term "uplift" potentials draw their origin from the first term of eq. ()2.9p 
while ip can be related, along the lines of [36], to an internal breathing mode. We shall keep 
this richer scenario in mind in the ensuing analysis, but in this paper we restrict our attention 
to the single-field setting. Taking into account the KKLT axion would complicate the analysis 
of cosmological perturbations, while as shown in ^ITj the early climbing behavior of the KKLT 
system is well captured by the scalar field alone. 

The single-scalar model affords interesting exact solutions for the simpler class of potentials 



whose analysis has a long history [211 [22l [23l [231 [25| [26] . To earlier results on general attractor 
trajectories and to the more complete analysis in [26], in particular in Russo's paper, [17J added 
more recently the recognition of some general implications of the climbing phenomenon, to which 
we now turn our attention. 

For the potential (12.10p . eq. (12. 8j) reduces to 




(2.9) 



V = 




(2.10) 




(2.11) 



that for < 7 < 1 admits the two distinct exact solutions (0 < r < 00) 



) 



(2.12) 



7 



and 



d(p 1 
~ch ~ 2 




7^ 



(2.13) 



Eq. (j2.12p describes a scalar field that emerges from the Big Bang climbing down the exponential 
potential, whereas the field if in (j2.13p initially climbs up, then reverts its motion and eventually 
climbs down. In both cases the scalar field eventually approaches in terms of the "parametric 
time" r the "limiting speed" 

^ (2.14) 



dr ^1-7 



2 



Notice that this limiting speed can be derived directly from eq. (12. lip , as for a simpler analogue 
of this system - a Newtonian particle in a viscous medium. It ceases to exist when 7 approaches 
one from below, and when this critical value is overcome the two classes of solutions in eqs. ()2.12p 
and (j2.13p leave way to a single class that describes a climbing scalar, with 



dip 1 
Ih ~ 2 



7+1 V2 ' ' / V 7 



(2.15) 



where now < r < , ^ [T7] . 



Considerations of this type capture important features of the two-term potential (j2.9p . whose 
dynamics is interestingly richer. If < 7 < 1, the asymptotic behavior for large positive ip is 
clearly dominated by the first term, which makes an initial climbing phase unavoidable, while 
the descent toward large negative values of 99 is dominated by the second term, so that ip is 
bound to approach asymptotically the limiting speed (|2.14p . Actually, when 93 comes close to 
terminating its ascent, the steep first term in eq. (12. 9p drives it to invert its motion somewhat 
abruptly, thus paying off less energy to cosmological friction than it would in the milder single- 
term potential (|2.10p . As a result, its descent proceeds readily at r-speeds larger than ()2.14p . 
which is thus approached from above. This qualitative feature will play an important role in 
Sections [3] and [H Finally, as discussed in [IT], if < 7 < the limiting speed (|2.14p translates 
into slow roll, so that even in the simplest manifestation of brane SUSY breaking, a one-field 
system with potential ()2.9p and 7 = the initial climbing is bound to inject an infiationary 
phase. 

An important exact solution of the single-exponential model, the Lucchin-Matarrese (LM) 
attractor [21], exists for < 7 < 1 and describes a scalar field that proceeds at the limiting 
speed (|2.14p . which is also the large-r behavior of generic solutions. In this case 

= V'o - -7==^ , -A = , (2.16) 

VI - 7^ Vl - 7 

where in the second expression we are omitting an additive constant that could be absorbed 
rescaling the spatial coordinates. In terms of the cosmological time t, defined by 

dt = e^d^ = , ^ dr (2.17) 

so that the metric (|2.2p takes the form 

ds^ = -dt^ + ei-^^^^fix-ffx , (2.18) 
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eqs. (I2.16P become 
1 



In 



7 



3 -i'^Mt 

2 7i 



7' 



A 



1 1 

7 7"^ 



In 



3 -f'^Mt 

2 7i 



1' 



and therefore the LM metric exhibits the power-hke behavior 



2 

37" 



dx- dx . 



(2.19) 



(2.20) 



2 713^2 ^ 

The LM attractor |2H [23l I25j thus describes a scalar field that actually slows down indefinitely 
in terms of the cosmological time t. This type of attractor solution plays a key role in the 
analysis of the CMB, and our task in the next sections will be to investigate how departures 
from it implied by the string-motivated picture of a scalar that climbs up the potential at early 
times can affect CMB observables. In a similar fashion, near the initial singularity, as the scalar 
climbs up the mild portion of the potential (j2.9p depending on 7, 

ip A — logr + ipo , (2.21) 

1 + 7 

1 

so that eq. (j2.17p implies that t ~ r i+t and 

,2 ,2 , 12 Mt 

ds^ -dr + 




, (l+7)(/3o 



3 1 + 7 



dx ■ dx . 



It is also useful to introduce the conformal time t], such that 

dt = e^-^^^^ drj , 



(2.22) 



(2.23) 



which is conventionally defined so that it approaches zero from negative values at late epochs, 
in terms of which the metric takes the form 



ds^ 

where the conformal scale factor is 

In particular, for the LM attractor 
?7 = 



'(r/) [-drj^ + dx- dx) , 



a(r/) 



e 3 



V6(l - 7^) 
Mil - 372) ^""^ 



37^ 



3^1 



1' 



so that 



ds- 




7 "Po 



drj + dx • dx 



(2.24) 
(2.25) 

(2.26) 



(2.27) 



We can now move on to correlation functions of gauge invariant perturbations around the 
background climbing solution of the potential (j2.9p . We begin with some basic aspects of cos- 
mological perturbation theory that are needed for our analysis, which are further elaborated 
upon in Appendix [Aj and then turn our attention to the so-called scalar and tensor spectra 
for the CMB anisotropies. These primordial spectra will form the basis of our comparison with 
observations in Section [SJ which we shall also reconsider in a WKB setting in Section O 
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3 Correlation functions off the attractor 



The models of the preceding section provide a natural mechanism to provide inflation its initial 
impetus. A phase of pre-inflationary climbing (beginning as fast roll) which subsequently injects 
a slow-roll epoch motivates the scenario wherein inflation commences even as the inflaton is off 
its eventual attractor. This section is devoted to exploring some general implications of this 
type of dynamics for the CMB power spectrum. 

We begin by considering scalar perturbations in terms of the Mukhanov-Sasaki (MS) variable 
v(r, ry), whose origin and applicability we briefly elaborate upon in Appendix [Al The starting 
point is provided again by the action (j2.ip . where the metric and the scalar field are both 
perturbed around some background solution {(j)o{r]) , a{r])) as 



and 

Defining the MS variable 

v(r,?7) 



(^(r,r?) = 0o(r/) + 5(/>(r,??) 

(l + 2'4){r,r])) dif + (l - 2 V'(r, r/)) dx • dx 



a{ri) 



a'(?7) 



(3.1) 
(3.2) 

(3.3) 



where a{r]) is the conformal scale factor of a metric in the form (|2.24p and "primes" denote 
derivatives with respect to the conformal time r], and letting 



z{ri) = o?{7]) 



a'iv) 



the quadratic action for v becomes 



5(2) 



d^x dr] 



(3.4) 



(3.5) 



The advantage of working with the MS variable v is that all reference to the original background 
is confined to the time-dependent mass-like term 



Wsiv) 



z"irj) 
z{ri) 



where z is defined in eq. ([3 



A similar action, with z(ry) replaced by 
zrp{ri) = a{r]) 
can be used to describe tensor perturbations. 



A 

e 3 



(3.6) 



(3.7) 



The MS variable is particularly convenient when comparing perturbations around two differ- 
ent homogeneous backgrounds that differ from each other perturbatively. For example, one can 
readily compare n-point correlation functions computed around the LM attractor of Section [21 
for which Ifl 

r?2 (1 - 372)2 ' ^•'•^^ 



In the LM case, this expression apphes for both scalar and tensor perturbations. 
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with corresponding correlation functions computed around a different background with 



= iy(0)(r?) + 6Ws{v) . (3.9) 

<^ 1 this can be done efficiently working in the interaction picture and 
treating 6Ws as a feeble interaction. The time evolution is determined by the operator 

U{'n,r]o) = Te''^^'^'^"' (3.10) 



In particular, if 



5Ws 



via the interaction Hamiltonian 



Hi = j d^k6Wsir]) |t?k(r?)|' , (3.11) 

with 

v{r,v) = [ ^k(r?) e^"^- , (3.12) 

J (27r)2 

where Wk(^) is the spatial Fourier transform of the MS quantum fielc^ v(r, r/) and all operators 
evolve, as in the standard Heisenberg picture, with the unperturbed Hamiltonian determined 
by eq. (j3.5p . To first order, if rj^ indicates the time at which the initial state was defined, the 
effect on the connected n-point correlation function is 

S{vi,Av) ■■■y^niv))in,in = -i / di]{[vi,^{r]) ...V]^„{r]), Hj{r])])in,in , (3.13) 

where the subscript on the expectation value indicates that, in the path-integral language, the 
expression ought to be evaluated resorting to the "in-in" Schwinger-Keldysh formalism |37j . 
Combining (13. 5p with given vacuum initial conditions and applying Wick's theorem one can 
compute for instance the correction to the 2-point correlation function as 



in, i 



z<53(ki+k2) r d?j6Ws{9j) 



(3.14) 



in terms of the Green functions for the LM attractor background of eq. (I2.27p . 



Gk{ri,?j) = -y^Hi'H-kv)Hi'H-kv) {v > v) , (3.15) 



where 



3 1 - 7^ 
2 1 - 37^ 



(3.16) 



The first-order correction to the 2-point correlation function picks up contributions from the 
integrand in eq. (13.141) that end up superposing features and oscillations if SWs{i]) has wide 
enough support in conformal time. It is instructive to see this in a concrete example, and to 
this end let us consider a step-wise perturbation. In order to simplify matters, let us also 



^In this paper quantum fields are denoted witli a "liat" in order to distinguisli from corresponding wavefunc- 
tions. 
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concentrate on the limiting case 7 = 0, which corresponds to standard de Sitter inflation, so 
that the corresponding Green function is simply 



If a constant 6Ws acts within the conformal time window (—??«, — ??/) with rjij > 0, the integral 
in eq. (j3.14p leads to an instructive deformation of the power spectrum, that as we discuss in 
the following section, can be defined in the present example through 

5\k + q)5PQ{k) = lim {r]Hf 5{0\vk{r])v^{r])\0) , (3.18) 



where 



5^{k + ci) Pgik) = k^ hm (0|Qk(r?)Qq(r?)|0) , (3.19) 



and Q{r,r]) is elaborated upon in Appendix |Al Therefore: 
6PQik) = kH^5Ws 



cos(2/c7?,) - -^r-^ - 7 cos(2%) ' ^ '^^ 



4 ' 2kr]i 4 ^ 'J' 2krif 



(3.20) 



If the actual Ws deviates from the attractor form of eq. (j3.8p for —rji < r] < —i]f, oscillations 
should show up in regions of the spectrum where \krji\ or \kr]f\ are of order one. In Section [5] we 
shall see this type of effects clearly in our numerical simulations. For the two-exponential system 
sizable oscillations draw their origin from early epochs close to the initial singularity, while 
additional features can be associated to a dip that, as we shall see in Section[6l develops in Ws{r]) 
around the endpoint of the climbing phase. Its depth reflects the entity of the bounce against 
the hard exponential, and indeed the oscillations that it brings about are more pronounced for 
scalar fields that penetrate more deeply into the steep region i;^ > 0. Finally, whenever \r]i\ is 
very large or \rif\ is very small the corresponding oscillations become essentially invisible, and 
indeed in Appendix |B] we shall present a class of exact power spectra with a sizable infrared 
suppression but no appreciable oscillations. 

Summarizing, one can anticipate the emergence of oscillatory features of scalar and tensor 
power spectra that typically present themselves even in the absence of initial particles. In 
the next sections we shall complement these arguments with numerical simulations and with a 
qualitative WKB analysis driven by analogies between the Mukhanov-Sasaki and Schrodinger 
problems. This last approach applies in principle even for large deviations from the attractor 
trajectory, and provides a simple rationale for the infrared suppression of scalar and tensor 
power spectra that is generic in models with an initial singularity. The overall IR modulation of 
the power spectra are also readily understood via spectral running arguments, which we discuss 
after presenting our results. 



4 The power spectrum 

In comparing the predictions of a particular model of the early universe to the observed CMB 
temperature anisotropics, one is required to compute the so-called dimensionless power spec- 
trum, defined as 

6'{k + ci)PQ{k) = k'{0\QUt)Qc,{t)\0) . (4.1) 
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In this section and the following one we shall work in cosmological time, and we shall also 
find it convenient to work with the rescaled MS variable Qk(^) of Appendix [X] that becomes 
independent of t shortly after horizon crossinj^. Qk{t) is a solution of the momentum-space 
counterpart of eq. (jA.lSp . 



fc2 Id 
+ Vs6 



.3i2 



■Xt) 



, 



(4.2) 



which describes perturbations around a background described by a{t) and (j)o{t). As in previous 
sections, "dots" indicate derivatives with respect to the cosmological time t, and V^^ is the 
second derivative of the potential with respect to (f) evaluated on the background solution, that 
here is denoted by 4>q. 

The 2~point correlation function is usually computed in the vacuum state for the background 
solution, for which ak|0) =0 for all k, where and Ok are creation and annihilation operators 
for the complete set of modes in terms of which the quantum field Q{t,:x) can be decomposed 



as 



Q(t,x) 



d^k 

(2^)3/2 

d^k 

(2^)3/2 L 



Ufc(t)e*''-^ak + ul{t)e 



-ik-x -^t 



(4.3) 



where the Uk are Klein-Gordon normalized solutions of (j4.2p . These modes will be different in 
general for different background trajectories ((/)o(t), a(t)), but computing (I4.ip in the vacuum 
state thus defined for our particular climbing trajectory yields the simple expression 



Pqik) = k^\uk{t)^' 



(4.4) 



which is to be evaluated after horizon crossing. For a Universe that has always been on the LM 
attractor, the potential term within the square bracket in eq. (14. 2p vanishes and one is left with 



Uk + 



1 



Uk + 



k^ 

_^ Uk 
t \ 37' 

to 



, 



(4.5) 



where we have chosen to so that a(to) = 1- The normalized modes corresponding to (|4.5p are 



Uk{t) 



3 7 ^ TT to ( t 



4(1 - 372) Vto 



3^Htl * 



to 



1 

37' 



V 



1 



37' 



(4.6) 



where hI^^ (x) denotes a Hankel function of the first kind of order u, which is defined in eq. (I3.16p . 
From this result and from the limiting behavior of the Hankel functions for small values of the 
argument [38j one can see that 



k-'\uk{t)\' 



k 



3-2u 



(4.7) 



for long wavelengths, so that the fiat spectrum of de Sitter space is recovered in the limit 7 — )• 0. 
On the other hand, the case 7 = corresponds to the spectrum displayed in fig. 1. 



'^"Furthermore (see also Appendix Qk(t) remains well defined throughout the climbing phase, and at late 
times can be converted into the comoving curvature perturbations that source temperature anisotropies. 
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Figure 1: Power spectrum for the LM attractor for 7 = COBE normahzed at k 

0.002Mpc~^ 



In evaluating the power spectrum for perturbations around a climbing solution, it is illumi- 
nating to allow for an initial state at Iq that may not correspond exactly to the vacuum specific 
to the climbing phase. For example, the modes describing the vacuum for a background that is 
exactly on the LM attractor are certainly different, as is the case for the Bunch-Davies vacuum. 
Therefore, one may be more interested in computing 

6^{k + q)PQ,^{k) = k^{^P\QUt)Q^{t)min,in ^ (4-8) 

where ip is an as yet unspecified state. For eigenstates of the number operator, or more generally 
tracing over a thermal density matrix, the result is a straightforward modification that reads 

PQ,Nik) = A:3|ufc(t)|2[l + 2nik)] , (4.9) 

where n{k) is the relevant occupation number in the k^^ mode. This follows directly from the 
fact that in such states (ak) = (Sj^) = while (akCip) is (5-function correlated^. More generally 
one can have an initial state that fits into neither of these two categories, but can still be 
parametrized assigning to the modes corresponding to k occupation numbers /3k- Specifically, 
on some fixed hypersurface, the initial data could correspond to all modes being in the vacuum of 
some background space-time other than what we are perturbing around (with associated mode 
functions Vk{t)). Such a state would appear to have finite occupation numbers with respect to 
the Q(t,x) quanta (with mode functions Uk{t)). On this hypersurface, the two sets of modes 
can be related via Bogoliubov transformations of the form 

Vk = cosh 6k Uk + e~^^^ sinh Ok ut , 
Vfc* = cosh 6k u\ + e * sinh 6k Uk ■ 

The initial state we are interested in is of the form I'i/;) = U{&)\0), where U is an operator 
that we shall specify shortly. The correlation function can therefore be recast in the form 

^'^These are special cases of the so-called "random phase" states, for which these properties also hold |39] . 
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in,in: SO that it Corresponds to performing a unitary transformation 
on the operator that has the effect of a Bogohubov transformation on the states. 

The quantum field producing particles of wavefunctions is 

Q'(t,x) = WQ{t,^)U = I [T4(t)e*'^-ak + ^^(t) e-*'^-^] (4.11) 



(27r)3/2 
(27r)3/2 



(4.12) 



where the two sets of creation and annihilation operators are related as 

6k = cosh 9k Ok + e*''* sinh 9^ alk ' 

6]^ = e"*^*-' sinh 9^ a_k + cosh 9^ ■ 

For transformations that only rotate modes up to some fixed momentum scale, the unitary 
transformation that formally effects this i j^ 

IJ{Q-^ = g-5 /^^^ [®fc^kS_k-©fcSk"-k] ^ (4-13) 

with 

Qk = Ok e'^'^ , (4.14) 
and one can choose 9^ > 0. The transformation U is defined so that 

U^a],U = 6k (4.15) 

and the result of starting with the initial data instead of is accounted by the unitary 
transformation U transforming the operator Q in g3]) as (5(t,x) ^ Q'(t,x) = U^Q{t,:x.)U. 

Relating the transformation (j4.10p to the standard Bogoliubov coefficients 

ak = cosh9k , Pk = 6"*'^'= sinh6'fc , (4.16) 

one can evaluate the 2-point correlation function in the state = U{Q)\Q) as 

(V'IQk(i)Qq(t)|V'>m,m = ( 0|Q'k(*) QqW |0>m,m = | Ffc | '(^^ (k + q) 

= |ufc|2{l + 2|/3fc|2 + 2 cos((5fc + 2Afe) \akPk\} + c^) , (4.17) 
where Uk = e^^''\uk\. In position space, the power spectrum is encoded in 

(V|Q^(t,x)|V^)i„,i„ = J ^^|nfe|2{l + 2|/3fc|2 + 2 cos(<5fc + 2Afc) \akPk\} , (4.18) 

so that 

PQ,^{k) = k^\uk{t)mi + 2\l3k\^ + 2 cos(5fc + 2Afc) \akpk\} ■ (4.19) 

As a result, the power spectrum picks up a modulation from the possible presence of initial 
particles. The effects of these initial excitations decay exponentially for all observable wave- 
lengths, but their possible presence should be accounted for when considering CMB observables 
that were generated within a short window of inflation. In Appendix [C] we show that for initial 
conditions that are natural for a Universe that emerged from a cold Big Bang dominated by 
the kinetic energy of moduli fields, the effects of any initial excitations are negligible. However, 
these effects may become significant if a large number of modes were in an initial state different 
from the one natural to the climbing phase. 



^^The Hilbert spaces that are related by the transformation in (|4.10[) are unitarily equivalent so long as Gfc has 
finite support, as will be the case in the examples of interest to us. If had support for arbitrarily large values 
of k {e.g. when is a constant), the two Hilbert spaces would be unitarily inequivalent [40j . 
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5 Numerics 




Fi gure 2: Background solutions for the potential ()2.9p with exponents 7 — and (/>o — —4 
(top), 7 = (middle two, with (f>o = 0,-4 respectively) and 7 = 0o = —4 (bottom). 
The dashed line in the third column indicates the LM attractor for the respective single term 
potentials V = ^e^^"^. 

In this section we present the results of our numerical integrations. They focus upon the 
simplest class of models involving a single scalar field minimally coupled to gravity that as we 
have seen in Section [2l combine an initial climbing phase with an eventual epoch of slow roll. 
The relevant potentials, 

Vicj)) = ^ (e^^ + e'^^A , (5.1) 



2 

were given in eq. ()2.9p . and we display them again here for the reader's convenience. They 
combine a critical exponential with another, milder exponential that is to sustain an eventual 
phase of slow-roll inflation, (p could represent the dilaton (or some other scalar field partly mixed 
with it), and therefore we shall be particularly interested in cases with a weak string coupling 
{(j) < 0), which are compatible with the climbing phenomenon. Starting with a particular value 
for (pQ, which is kept as a free parameter, the initial speed for (p at very early times is determined 
by (I2.13|) . to be expressed in cosmological time t, for any given slow-roll exponent 7. Lettinj^ 
= we integrated the background equation of motion for (p and the results are presented 



Where we remind again the reader that in restoring k in eq. (|5.1|) . the dimensionful pre-factor of the potential 



would read jj^, so that M has indeed the dimension of a mass once we set k — I. 
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in fig. El Our integrations begin at tM = 0.01, where the scale factor is normalized to a = 1, 
with a — )• as t — )• representing the initial singularity. 

Given that H^t^ — )• on the attractor, the third column of fig. [2] shows that the solutions 
approach the attractor asymptotically long after the effective onset of inflation, as can be seen 
from the second column that displays the slow-roll parameter]^ 

H ^ -§2^ (5-2) 
- (5.3) 

where (pQ is the background solution and V^^ is the second derivative of V with respect to (p. 
Note that the accelerated expansion (e<^ < 1) begins within one "string" time tM ~ 1 in all 
examples, not much after the climbing phase has ended (left column of fig. [2]). Although r]^ takes 
slightly longer to converge to its attractor value, its relative size does not preclude the effective 
onset of accelerated expansiorj^. We note in passing the brief spurt of de Sitter inflation (e^ = 0) 
around the inflection point. Although the solution closes in on the attractor only at rather late 
times ~ 0(10^), inflation off the attractor begins almost immediately. Although inflating 
off the attractor is far less effective (in the sense that far fewer e-folds of expansion occur in a 
given interval of time off the attractor than on it), a significant number of e-folds can still be 
attained within a relatively short time (as measured in units of M). 

We have computed the power spectra for a range of slow-roll exponents in figs. [3] and HI 
We note that, in addition to superposed oscillatory features, at long wavelengths there is a 
marked suppression of power in the scalar spectrum and a corresponding, if less pronounced, 
enhancement of power in the tensor spectrum. We have adapted these results to the best fit 
spectra generated by equivalent solutions that are perpetually on the attractor, determined by 
the criterion that for the shorter wavelengths (spanning the majority of the scales that we can 
access with current observations) the two types of spectra coincide 

In each case, the plots represent comoving wavelengths that have exited the horizon over 
approximately 12 e-folds of inflation. Given that the range we access in WMAP data constitutes 
a mere window of 6 e-folds, we see that we are at most privy to a subset of the modes represented 
in these plots. Although we would have to be exceedingly privileged to access wavelengths that 
exited right as inflation begar0, we see that one can still expect to flnd some trace of the 
climbing phase if the largest modes that are accessible to us left the horizon no later than 6-7 e- 
folds after the climbing phase. This is due to the fact that although inflation effectively started 
right away (within one string time tM ~ 1 of the end of the climbing phase), the eventual 
approach to the attractor takes far longer. This is illustrated by the third column of fig. [2l 



An alternative parametrization for the second slow roll parameter is sometimes taken to be r;^ = which 



is related to 77^1, above during slow roll as r;^ = £</> + ??,/). We prefer the definition (|5.3|l . as this parameter is well 
defined throughout the range of our integrations. 

'^^We remind the reader that 77^ is only required to be small if its value is taken to be approximately constant 
during inflation. One can readily accommodate temporarily large values of 770 consistent with persisting inflation 
provided it becomes negligible shortly thereafter, since the existence of an accelerated phase only requires that 
be small at any given time. 

'^^Short wavelengths in this case are defined demanding that the climbing phase be sufficiently outside our 
present horizon and that at some convenient pivot scale (e.g. 0.0805Mpc~^) as defined by the horizon, the two 
spectra be COBE normalized. For the purpose of our plots, we have taken this scale to define the comoving 
wavelength k = 10*M. 

^'^As we have stressed in the Introduction, however, this possibility is to be rightfully considered [11] in the 
event that it may explain some large-scale features that do not fit comfortably within ACDM models. 
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Figure 3: Scalar (left) and tensor (right) spectra for 7 = ^^I'/'o = ~4 (top), and 7 = ^^I'/'o = 
—4 (bottom)- best fit at short wavelengths with attractor spectra for ns = 0.944, rij- = 0.9375 
and ns = 0.995, ut = 0.97, respectively. The comoving wavenumber k is in units of M. The LM 
attractor would have led to ns = nr = 0.9375 for 7 = and ns = nx = 0.972 for 7 = -y= . 



We have also evaluated the multipole moments of the angular power spectrum for the climbing 
solutions corresponding to 7 = and respectively, and we have compared them to the 
corresponding best-fit power-law spectra of figs. [3] and [H For this, we simply used the fact that 
at the largest angular scales {i < 30), to an excellent approximation [5T] 



Vr 



(5.4) 



where ji is a spherical Bessel function and ijq — r]r denotes our present comoving distance to last 
scattering surface. In this fashion we find again that if inflation had started within 6-7 e-folds 
of our present horizon exit, the climbing phenomenon would bring about a noticeable drop in 
power at the largest angular scales that would become more significant the closer the climbing 
phase were to the exit of our current horizon. 

We note that the drop in power for a range of parameters (depending on the multipole in 
question) becomes more significant than the uncertainty implied by cosmic variance, 



2^ + 1 ' 



(5.5) 
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Figure 4: Scalar (left) and tensor (right) spectra for 7 = with (pQ = (top), and 0o = —4 
(bottom) - best fit with attractor spectra for ns = 0.979, rix = 0.958 and ns = 0.975, n-r = 0.958 
respectively. The comoving wavenumber k is in units of M. The LM attractor would have led 
to ns = riT = 0.96. 



for each of the multipoles represented in fig.[5l if A; = Hq had exited the horizon during inflation 
no more than 2 e-folds after inflation began. Therefore a statistically significant drop in the 
power at long wavelengths can result if the climbing phase occurred within 3 1/2 e-foldJ^ of 
the horizon exit for the largest scales that are presently accessible to us. It is tempting to infer 
that this may be at least partly responsible for the apparent large-angle suppression of power 
in the CMB [16j. 

It is informative to compare the best-fit values of ns and nx obtained in figs. [3] and H] 
(indicated in the respective captions) to those of corresponding solutions that would have always 
been on the LM attractor corresponding to the particular exponent 7. Up to the expected 
higher-order corrections coming from the time variation of e^f, and r]^, they are in reasonable 
agreement at short wavelengths. Furthermore, one can readily understand the gross features of 
scalar and tensor spectra starting from the spectral relations for solutions that remain on the 

^*In all examples that we have displayed, with 7 = and about 1.5 e-folds of expansion resulted before 
slow-roll was effectively attained. 
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Figure 5: Changes in the scalar and tensor multipole coefficients for 7 = with respect 
to featureless spectra with ng = 0.944, nr = 0.9375 (top), and for 7 = with respect to 
featureless spectra with ns = 0.995, tt-t = 0.97 (bottom). 

attractor [42] : 

ns - ^ = 2{r]^ - 3e^) , (5.6) 
UT - I = -2e^. (5.7) 

For the single-exponential potential (12.101) . on the attractor, = 87^, ry^ = 67^ and so — 1 
= TT-T — 1 = — 67^. Reconsidering now eqs. ()5.6p and ()5.7p while allowing for spectral running (i.e. 
allowing for ns = ns{k),nT = nT{k)), and deriving these from the values obtained for e(f,{tk) 
and r)^{tk), evaluated at the time at which comoving wavenumber k exits the horizoil^. one can 
arrive at the running spectral indices plotted in fig. [6l We note that in each of the plots, the 
spectral indices run towards their attractor values for short wavelengths (as the solution settles 
towards the LM attractor), but deviate considerably for longer wavelengths. We can now readily 
understand the long wavelength modifications of scalar and tensor spectra - from (j5.7p we see 
that as monotonically decreases toward its asymptotic attractor value as it emerges from the 
climbing phase, n-r — 1 approaches —67^ from below. This means that for longer wavelengths, 
the spectrum had an even redder tilt, hence accounting for the convexity of the envelope for 
the modulations of the tensor spectra in figs. [3] and HI consistently with the enhancement of 
power at the longest wavelengths. For the scalar spectra, from eq. (|5.6p we find that the slower 
convergence of ry^ to its attractor value compared to results in a positive (blue) spectral tilt 

^^Given the fact that — f]^ = 0{e^, ri^,rj^€^), allowing for such a spectral running will be accurate up to 
errors quadratic in the slow-roll parameters. 
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Figure 6: Running spectral indices for 7 = (top) and 7 = (bottom). The comoving 
wavenumber k is in units of M. 



for longer wavelengths, which eventuahy asymptotes to its negative (red) attractor valuq^- This 
accounts for the concavity of the low-frequency portion of the scalar power spectrum, and hence 
for the wide suppression of scalar power at long wavelengths. 

One can also understand the qualitative effects of starting inflation off the attractor, but now 
initially descending, which would be possible if the first exponential in eq. (jS.ip were replaced 
by a slightly milder one. In this case the spectral index ns{k) would start out again larger than 
one, but cosmological friction would make it converge to one more rapidly, resulting in a less 
relevant suppression of the scalar spectrum that would be confined to much longer wavelengths 
in comparison with the initially climbing solution. All in all, the effects of pre-inflationary 
dynamics are far more significant for a solution that initially climbs, since the scalar field comes 
close to the attractor trajectory much later in this case. 

Although we find a noticeable suppression for the lowest multipoles if the largest scale we 
observe today exited the horizon within 6-7 e-folds of the climbing phase, this drop would 
be statistically significant [^3 only if horizon exit occurred within 3-4 e-folds. However, the 
consistency of the suppression over a wide range of multipoles ought to boost its significance 
even if horizon exit occurred at later times. Furthermore, if for various parameters of the 
potential and initial conditions for the inflaton, the superposed oscillations persisted to shorter 
wavelengths to a significant enough degree, the combination of the two effects might be easier to 
discriminate and constrain in the data. We leave to a future study a full analysis of any improved 
fit to WMAP data that may be offered by a climbing phase. However, in merely setting out 



^"Note the coincidence of the inflection points of the scalar spectra in fig. |3]with the transition of the spectral 
index from blue to red in flg. [G] 

^'^In the sense that the drop in power for the low multipoles would surpass the threshold obscured by cosmic 



variance. 
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to understand the generic consequences of a class of pre-inflationary dynamics inherent in an 
interesting class of string constructions, we find it encouraging that large-angle anomalies in the 
CMB [16] may find a ready and natural explanation. 



6 The Schrodinger analogy 



Further qualitative insight into the modifications of the CMB power spectrum can be attained in 
general via WKB methods. The MS equation that results from (j3.5p has in fact the Schrodinger- 
like form ^ 

(4^ + pHk,7])] vkin) = , (6.1) 



where 

„2/L „^ _ ;,2 



p'{k,ri) = e - Ws{r]) (6.2) 
and where the MS "potential" Ws is defined in eq. ([35]) El- 



Notice that some key information is generally available on Ws for this type of systems, which 
determines completely its asymptotic behavior at the initial singularity and at late times. One 
can indeed show, starting from the results in Section [21 that 

Ws - T 7 ^—TT7 , (6.3) 

where — r^o (with r/o > 0) corresponds to the initial singularity and is determined by the condition 
that the conformal time ry tend to zero for large values of t or r. In addition the systems of 
interest, with < 7 < approach an eventual slow-roll phase, so that 

, (6.4) 

where v is defined in eq. ()3.16p . since Ws is to attain the limiting LM behavior (j3.8p . These 
two limiting behaviors make it inevitable for Ws to cross the real axis for a certain value of the 
conformal time in the interval {—r]o,0), as can be seen in fig. [71 As we shall see, this generic 
structure suffices to bring along a generic depression of the power spectrum of CMB fluctuations 
in the infra-red in models with an initial singularity. 

There is therefore a striking analogy with one-dimensional Quantum Mechanics, but for 
an important proviso. Systems describing an eventual inflationary phase map into an initial- 
value Schrodinger-like problem featuring an impenetrable barrier that, as such, overcomes any 
available "energy" k ^ for small negative values of rj, as can be seen in fig.[7l In ordinary Quantum 
Mechanics, a barrier of this sort would result in total reflection, and thus in a stationary wave 
in the classically allowed region, accompanied by a decaying mode in the classically forbidden 
region. In the Mukhanov-Sasaki problem, however, in the intermediate region where Ws{ri) is 
small, Vk{r]) and its complex conjugate ought to be close to incoming and outgoing plane waves 
if close to the initial singularity the system is to lie in a Bunch-Davies-like vacuum. In this 
fashion, beyond the barrier these types of solutions leave way in general to combinations of a 
decaying and a growing mode, the latter of which readily dominates. In a cosmological setting. 



^^As we have seen, the corresponding MS equation for tensor perturbation is identical, up to the replacement 
of Ws with the "potential" Wt determined by eq. (|3.7p . 
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Figure 7: Ws{ri) for a climbing solution in the potential (j2.9p with 7 = 0.1 (continuous red), 
and W^'^\ri) for the corresponding LM attractor in the potential (I2.10p (dashed blue). 



this very behavior reflects itself in the growth of perturbations after horizon crossing. The gross 
features of the phenomenon are well captured by the conventional WKB matching relations, 
which in the classically forbidden region lead to 



t'fc(-e) 



1 



exp 



\p{y)\ dy 



(6.5) 



where the small positive number e << 1 identifies the end of inflation, and where at early enough 
initial times an oscillatory behavior like 



1 



- ikrj 



/2I 



(6.6) 



would correspond to the Bunch-Davies vacuum of de Sitter inflation. In eq. (|6.5|) — ??* is the 
classical turning point of fig. [3 such that Wsi—f]*) = k'^, and one is interested in the leading 
behavior for e — )• 0, after several e-folds of inflation. 

For a scalar field on the LM attractor corresponding to the milder exponential potential of 
the previous section, 

(0.7) 



V 



with 7 < eq. (16. Ij) reduces to 



where v is defined in eq. (13.16P . One can appreciate the nature of the leading WKB approxima- 
tion comparing it, for the LM attractor, to the exact solutions of eq. (j6.8p . As is well known, 
these can be expressed in terms of Hankel functions as 



Vkiv) = e^^^^'^^^) H^J\-k^) 



(6.9) 
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so that the hmiting behavior for rj ^ —e [38l 




1 



— V 




(6.10) 



results in the power spectrum 



~ k 



3 2 



k 



(6.11) 



z 



with the characteristic 3 — exponent. 

The essential features of this expression are qualitatively, albeit not exactly, captured by the 
leading WKB formula ()6.5p . which up to overall phase factors yields 



The WKB amplification factor of eq. (I6.12P thus recovers the qualitative features of the power 
spectrum of eq. ()6.1ip . but getting the precise slope would entail the replacement of v with D, 
taking into account higher-order corrections El- 
Together with the definition (|6.11l) of the power spectrum, eq. (16. 5p readily accounts for a 
generic suppression of the low-frequency tails of scalar and tensor power spectra in models with 
an initial singularity. For any given k, if one begins in the generalization of the Bunch-Davies 
vacuum for negative 77's in the wide flat region of Wsif]), the amplification of Wfc(r/) at horizon 
crossing reflects the area below |j3(?7)|, or equivalently the contribution of regions of conformal 
time where Ws{i]) > k"^ (see fig. [7]). As we have stressed, Ws is bound to be negative at early 
times, so that the amplification factor cannot grow beyond limit for small values of /c, and 
precisely in this region the overall k^ of P-jzik) takes over. 

On the other hand, the WKB formula (j6.5p does not capture an interesting feature related 
to the dip present in fig. [71 which belongs to a region where Ws is negative and is roughly 
centered around the conformal time when the scalar inverts its motion. One is thus led to expect 
that oscillations be present in the power spectrum, as for the transmission past an attractive 
well in Quantum Mechanics. This reasoning is along the lines, and actually equivalent to, the 
perturbative argument presented in Section [3l As we saw in Section [5l oscillations are indeed 
generically present in the power spectra. 

Fig. [7] displays typical forms of Ws{r]) for the two cases of interest, an LM attractor and the 
climbing solution for the two-exponential potential of eq. (j2.9p that eventually approaches it. 
The climbing solution emerges from an initial singularity at a finite negative value of r], where 
Ws{ri) has a vertical asymptote, to then approach for — t- the curve W^^\rj) of eq. (|3.8p . 
Actually, Ws{i]) for the climbing scalar lies always below W^^\rj) and comes close to it rather 
late, so that for relatively small values of k the integral in eq. (j6.5p is reduced with respect to its 

■^^The apparent discrepancy indeed disappears if one takes into account higher-order WKB corrections, which 
recover the expansion of u, expressed as J-y^l + j^, in inverse powers of T^, as the reader can verify for the next 
few orders. The successive orders of the WKB expansion are in fact sized by the relative variation of the potential 
within spatial distances comparable to a wavelength, which is smaller for larger values of u or u. 





where 




(6.13) 
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Figure 8: H^t(??) for a climbing solution in the potential (j2.9p with 7 = 0.1 (continuous red), 
and W^^\ri) for the corresponding LM attractor in the potential (I2.10p (dashed blue). 



attr actor value for two conspiring reasons. First, a lower Ws{ri) for the climbing scalar implies 
that the lower end —r]* of the integral in ()6.5p lies to the right of the corresponding value for 
the attractor, and in addition the area below the curve, and thus below |p(r/)|, is also reduced. 
All in all, one can thus see by inspection that the WKB amplification factor in eq. (I6.5p is 
reduced for the climbing scalar, and more so for smaller values of k. On the other hand, for the 
small negative values of rj that are relevant for large values of k, Ws{r]) approaches from below 
W^^\ri), and therefore the power spectra of climbing scalar and attractor solutions ought to 
merge in the ultraviolet region. Conversely fig. [8] indicates that, tracing it backwards from the 
high-frequency region where it falls on the LM form of eq. (j6.1ip . the power spectrum of tensor 
perturbation ought to be initially slightly enhanced in the infrared region with respect to the 
attractor solution, before being finally suppressed at very low frequencies as Wt(??) crosses the 
real axis. This is all along the lines of Section [H and actually the relative behavior of scalar and 
tensor perturbations can be foreseen rather clearly via the WKB method, comparing the two 
expressions 



2 



dr 



1 + 



dip\ 
dr I 



V 

=. + 9 



17 ( f )' 



1+ 



and 



Wt ^ = —63-^ 

a 3 



(6.14) 
(6.15) 



which in their turn determine the corresponding spectra via the WKB amplifications factors. 
Here a is the conformal factor of eq. (I2.25P and "primes" denote as usual derivatives with respect 
to the conformal time ij of eq. (j2.23p and one is interested in the final descent, where the second 
exponential in eq. IK9\\ dominates so that 27 F and K^<^ y. 

One can verify that the ratio which determines via eq. (j6.5p the corresponding effects 
on the power spectra, equals one for the limiting speed (j2.14p and is smaller for the larger 
values that, as we have argued in Section [2l characterize the descent in the two-term potential. 
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Actually, for ^ close to the attractor value (|2.14p one can show that 

Wt (2-372) ^^dT ^/T^) ^ ' 

and therefore the effect presents itself symmetrically about the attractor speed and is in principle 
sizable for the relatively small values 7 < that grant an eventual slow roll. 

In the numerical investigations described in Section [5l scalar and tensor perturbations were 
nicely behaving along these lines. There we had also reconsidered our findings in terms of 
spectral running as the inflaton comes close to the attractor, and the arguments of this section 
provide complementary information to the same effect. In Appendix [B] we shall present some 
analytic models for Ws whose exact solutions capture the qualitative features displayed by the 
power spectra of figs. [3]and[ll aside from their oscillatory patterns. 



7 Conclusions 

In this paper we have analyzed the imprints on the CMB of a non-standard phase of cosmological 
evolution that manifests itself before the onset of slow-roll inflation in a class of string models 
where SUSY is inevitably broken at the string scale. In this phase, the only consistent solution 
for the scalar field which connects to an initial curvature singularity (the cold Big Bang), forces 
it to be initially climbing up an exponential potential [TTj. In the simplest two-exponential 
potential (j2.9p . this scalar field quickly reverses its motion and begins to slow-roll inflate, even 
as it remains rather displaced off its attractor solution in held space. This scenario was shown in 
p!7] to be inevitable for steep exponential potentials, whose logarithmic slope reaches (in absolute 
value) a critical limit that, interestingly, is attained both in the simplest ten-dimensional string 
model with "brane SUSY breaking" [29] and in the four dimensional KKLT scenario [31] with 
an F-term uplift [l3]. In models where a steep potential term appears in combination with 
another slow-roll exponential, the climbing phenomenon is followed by a subsequent slow-roll 
inflationary regime. Whereas we discussed the climbing phenomenon within a specific context, 
these results refiect a more general Kasner-like behavior that can be exhibited working directly 
in cosmological time close to the initial singularity, proceeding along the lines of the BKL 
analysis in General Relativity [44] . In this more general setting, one can see that the climbing 
phenomenon continues to occur in more general scalar potentials containing steep portions that 
constrain the behavior right after the singularity. Higher-curvature string corrections might also 
be compatible with the climbing phenomenon, at least in special classes of constructions, but 
we do not have conclusive arguments to this effect at the present time. 

We have shown that if such a climbing phase had occurred within six to seven e-folds to when 
the largest scales accessible to current observations exited the horizon, specific imprints mani- 
fest in the CMB power spectra. In addition to long-wavelength suppression and enhancement 
of the scalar and tensor primordial spectra respectively, both are generically imprinted with 
superposed oscillations. We provided some semi-analytic arguments for these features based 
on the in-in formalism for computing cosmological correlation functions, on a WKB analysis of 
fluctuation modes and on special analytic models for the MS potentials Ws and Wt- Further 
indications were provided by the spectral running implied by inflating off the attractor. This 
modiflcation of power at long wavelengths translates directly into a reduction in low multipole 
moments of the scalar angular power spectrum, which could become more significant than cosmic 
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variance limits if the climbing phase occurred within three to four e-folds of when the largest 
scales accessible to us exited the horizon, and in a corresponding but milder increase in low 
multipoles of the tensor power spectrum. It is of course difficult to argue that the suppression of 
the scalar power spectrum at large wavelengths is an unambiguous predictions of the climbing 
phenomenon, but for instance for descending solutions, an option that appears disfavored in 
String Theory but would nonetheless be available if the first exponential in (j2.9p were replaced 
by another lying below the critical slope, or in the absence of the "hard" exponential brought 
about by String Theory, the consequences for CMB observables would be far less dramatic. In 
order to obtain inflation successfully in the class of models we consider, the mass parameter M 
defined in Section [5] should be of order 10^^ — 10^^ GeV, which is reasonably achieved in the 
SUSY breaking mechanism of |29[ [20] that is closely linked to the string scale. In a 4d KKLT 
context, this scale would reflect an energy scale of the uplift, or an appropriate combination 
of string scale and vev of the axion-dilaton fields in a compactification of the ten-dimensional 
"brane SUSY breaking" setup. 

It is tempting to ascribe at least part of the large-angle anomalies in the CMB [16], and 
in particular the evident lack of scalar power at the largest observable wavelengths [l5], to a 
pre-inflationary dynamics of this type. We hope to have convinced the reader that the class of 
models that we have analyzed, well motivated by String Theory, can provide an aesthetically 
compelling scenario for the Early Universe that not only gives a natural impetus for inflation 
to begin, but can also offer potentially observable predictions for the CMB that might allow a 
better fit to existing data. 
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A COSMOLOGICAL PERTURBATIONS 



Let us briefly review the basics of cosmological perturbation theory that are required for our 
investigations. Our starting point is the action for gravity coupled to a canonically normalized 
scalar field 



S 



d —g R + d X\/—g 



and we consider perturbations around the background solution for the scalar field 

(/>(r,7?) = (j)o{r]) + 5(t){r,r]) 

and for the line element 



ds^ 



(l + 2'i/'(r,r/)) + (l - 2V'(r,r/)) dx' dx' 



(A.l) 



(A.2) 



(A.3) 



We work around an arbitrary background solution (j)o{i]), a{rj) in the so-called longitudinal gauge 
|46j . where the only independent metric perturbation is ^/'(r, ry), with r/ the background conformal 
time. Substituting this ansatz into the action and making the crafty field redefinition 



v(r,?7) = a{ri) 



a'(j]) 



where primes denote derivatives with respect to conformal time, and letting 
the second order variation of the action becomes simply 



5(2) 



- / d-'x drj 



(v') 



/\2 



(A.4) 



(A.5) 



(A.6) 



We note that eq. ()A.6P is formally the action for a scalar field on a Minkowski background 
with a time dependent mass-like term. v(r, rj) is commonly known as the Mukhanov-Sasaki 
(MS) variable, and is related to the comoving curvature perturbation TZ{r, rf) according to 

v(r,r?) = z(77)7^(r,r?) . (A.7) 

Eq. ()A.6P is valid for any consistent conformally flat solution of the coupled system of eq. (jA.ip . 
but it is not the only relevant form of the action for the perturbations since one could have 
chosen a field redefinition different from ()A.4p . Thus, working directly with IZ and switching to 
cosmological time, eq. ()A.6P can be recast in the form 

= -\ j d'x^g[^)"d,nd^n, (A.8) 

which would be the standard action for a massless scalar field in the given background spacetime 
sourced by the dynamics of (j)Q, were it not for the time-dependent normalization of the kinetic 
term. For our purposes, another convenient choice is the rescaled MS variable 

Q(r,t) = , (A.9) 
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in terms of which the action can be cast in the form 



5(2) 



^ J d^x^d^Qdf'Q - I j d^x^g [V^^ + Q{t)]Q^ , (A.IO) 



with 



nit) 



1 d /a^<j)l 



dt\ H 



d f2V 



(A.ll) 



where the equivalence between the last two expressions follows from the background equations 
of motion (|2.6p and (j2.7p . Alternatively, one can recast the action in a form involving the two 
(not necessarily small) parameters, 



dr J 



1+ 



2 ' 



V 



(A.12) 
(A.13) 



The end result. 



5(2) 



d^x 



-g { d^Q d^'Q + 



H 



+ (6e^- 2e| + 3/?^ 



Q' 



(A.14) 

is valid for arbitrary backgrounds. The special case where e,^ ^ 1 implies that the background is 
instantaneously undergoing slow-roll inflation, whereas the smaller r;^ is, the longer this period 
of slow roll can be taken to last. Typically, a long inflation (defined as more than 60 e-folds 
worth) requires r/,^ ^ 1 as well. The advantage of working with the rescaled MS variable Q{v, t) 
is that it is constant on super-horizon scales, and remains well defined throughout the evolution 
of (po, including during its early climbing phase. 

Finally, the resulting equation of motion for the Fourier modes Qk(0 of the operator Q{r,t) 
of Section [H takes the relatively simple form 



+ "J-" — 77' + 



dt 



+ Vm + 



d_ 

It 



2V\ 



Qk = , 



(A.15) 



which upon use of the background equations of motion, results in the equivalent form ()4.2p . We 
make use of the action in terms of both v and Q perturbations in the main body of the paper. 



B Analytic models for Ws and Wt 

As we have stressed in Section [6l there are interesting analogies between the one-dimensional 
Schrodinger equation and the MS equation in the form (|6.ip . provided one identifies /c^ and Wg 
with energy and potential for the quantum-mechanical problem, together with a key difference, 
since in the latter case one is solving an initial-value problem. This has the key consequence 
that the infinite barrier resulting from the eventual approach to the LM attractor crystallizes 
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-0 



Figure 9: Comparison between attractor Ws (dashed blue), two-exponential Ws (continuous 
red) and the approximate form of eq. (2) (dotted brown). 



and largely amplifies initial fluctuations. Elaborating on the analogies between the two settings 
can provide a better qualitative grasp of the power spectra associated to various cosmological 
settings. In this Appendix we investigate some relatively simple analytical forms for the poten- 
tials Ws and Wt that suffice to capture the gross features of the power spectra in figs. [3]and|H 



Leaving aside momentarily the Big Bang singularity, there is an interesting class of deforma- 
tions of the attractor MS potential 



1 



VF(0)(ry) = (B.l) 

where 



jj2 



2 1 _ 2-37 
" "4 - (1-37^)2 ' (^-'^ 



that captures the gross features of the actual MS scalar and tensor potentials for the two- 
exponential system over most of the cosmological evolution. Referring for definiteness to Ws, 
the relevant potentials read 



Ws- ^ 



^2 



(B.3) 



and combine the proper LM late-time behavior, a single zero and an almost flat region. These 
potentials depend effectively on a single interesting parameter, c, since —rjo is simply the location 
of the zero, and have the virtue of preserving the attractor behavior of eq. (16. 4p as 77 — )• 0~. 
The relevant range is c > 1, in order that the potential tend to zero from negative values as 
r] —7- —00, which for one matter simulates the almost flat negative plateau present in the actual 
forms of Ws and Wt- An important transition clearly takes place at c = 2, since for 1 < c < 2 
the potential of eq. ()B.3p lies below the LM potential at late times while the opposite is true 
for c > 2. For the sake of comparison, fig. [9] displays the actual Ws determined by the two- 
exponential potential (12. 9j) with 7 = 0.1 and ipo = 0.5, the corresponding attractor Ws and a 
good approximation provided by eq. ()B.3p . obtained from eq. ()B.3P with c = 1.1285 and r/o = 1. 
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For the potentials (1B.3P the MS equation reads 



(f Vk 

or, letting 



+ 



:i - c) u 



2 _ 1 
4 



) (2-c)(.^-i) 



Vk 



{krjor + (c-l)(z.2--) , 



+ 



^0 



(2 - c) (..2 - 1) 



■2 _ 1 
4 



. 



The additional substitution 



makes p positive in the relevant range and turns eq. ()B.6P into 



d"^ Vk 
dp^ 



+ 



Finally, letting also 



1 + 



1 



(2 - c) - i) 



Ap 



Vk 



(B.4) 

(B.5) 
(B.6) 

(B.7) 
(B.8) 



1 



L{L + l) = iy--^L = u- -, a 



(2 - c) - 1) 



2A 



(B.9) 



eq. (IB.Sp takes the standard form of the equation for Coulomb wave functions discussed for 
instance in Chapter 14 of [38], 



d^Vk 
dp^ 



+ 



2 a L{L + l) 



P 



Vk 



. 



(B.IO) 



Taking into account the unit normalization for the Wronskian in terms of r], the relevant solution 
for us is the combination (see eq. ()B.6p . where in the flat region A^/?7q plays the same role played 
by k"^ in the MS equation for the attractor) 



Vk 



(Fl + iGL) , 



(B.ll) 



where Fl is the regular Coulomb wavefunction and Gl the irregular one. Near p = is 
dominated by the second term, so that the k dependence of the power spectrum is captured by 
the limiting behavior of 

k^ 



Pnik) ~ ^ 



(B.12) 



as e tends to zero. This can be extracted from eq. (14.6.7) of [38j and from the power-like 
behavior of the modified Hankel function K2L+1 near the origin, and reads 



Pn{k) 



^3 gTra 



\T{L + l + ia)\^ A2i+i 



(B.13) 
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Figure 10: Analytic scalar (red) and tensor (blue) spectra vs attractor spectrum (dotted). 



or, in terms of our original variables, 

(k r]o)^ exp 

Pnik) ~ 



-(f-i)('-^-i) 



. (B.14) 



This expression combines the typical attractor behavior P{k) ~ k^'"^^ at large frequencies 
with modifications at long wavelengths that, for any given 7, depend both on t/Qj which sets the 
scale of the phenomenon, and on c. The scalar power spectrum for the preferred potential Wg 
corresponding to fig. [9] is the red curve displayed in fig. [TOt while the blue curve is the tensor 
powers spectrum corresponding to c = 3 and ryo = 3 that results from a similar fit for Wt of 
fig. El The infra-red depression present in both spectra and relative enhancement that follows it 
in the latter curve should be compared with the gross features of the scalar and tensor spectra 
of figs. [3] and [H One can also verify that the leading WKB approximation (16. 5p captures the 
behavior of eq. ()B.14p for large v. 

Relatively simple analytic potentials capturing also the Big Bang singularity can be obtained 
at the price of turning to the less familiar and more complicated Heun functions |47j . Suffice 
it to say, here, that just like the Bessel and Coulomb cases describe systems with the late-time 
regular singular point [r] = 0) and one irregular singular point at infinity, the relevant confluent 
Heun equation possesses two regular singular points (corresponding to the initial singularity and 
?7 = 0) and an irregular singular point at 00. 



C Initial state effects 



In this Appendix we would like to elaborate on the effects of initial particles on the power 
spectrum derived in eq. (|4.19p : 

Dk = {l + 2\pk? + 2 cos(5fc + 2Afc) |afc/3fc|}. (C.l) 

Clearly, for isolated resonances of specific momenta, {i.e. for an incoherent superposition of 
a finite number of single-particle states), the power spectrum would acquire small localized 
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"bumps" . As these particles are rapidly diluted by the expansion of the Universe, their contri- 
bution to the power spectrum over an observable range of wavelengths would decay exponentially. 
In considering their presence, one should therefore justify or derive their appearance in the first 
place. Generically, one would arrive at Bogoliubov coefficients that would be too negligible to 
impart much more than a small amount of "noise" over the power spectrunJ^. 




10 100 1000 10* 10' W 10 100 1000 10* 10' w 



Figure 11: Modified power spectra for 7 = (left), and 7 = (right). Dotted lines are the 
climbing scalar power spectra, dashed lines are the corresponding attractor spectra. 

In the scenario that we considered in the main body of the paper, the Universe begins with a 
"cold Big Bang" , wherein it emerges from an initial singularity dominated by the dynamics of its 
lightest relevant degrees of freedom ~ scalar fields representing various moduli in String Theory 
evolving in the potentials that have been generated for them. The fluctuations of the moduli 
around their (rapidly) rolling solutions can naturally be taken to be in their vacuum state, with 
the eventual hot Big Bang arising at the end of inflation, when the Universe reheated via the 
production and subsequent thermalization of Standard Model particles from the energy stored 
in the inflaton. However, one might ask what effects the CMB observables would display if the 
fluctuations around the climbing background had been in the vacuum state natural to some 
other solution. Specifically, what if the Universe started out in a radiation dominated epoch, 
only for a climbing scalar to condense and take over the evolution? Although this is admittedly a 
somewhat concocted situation, it is illustrative to consider its consequences, since the initial state 
induced on the climbing phase will have non-vanishing occupation numbers for long-wavelength 
modes that one can calculate. In general, if at an initial time Iq the perturbations were in the 
vacuum state relevant to a Universe expanding as , computing the Bogoliubov coefficients 

in eq. (|4.19p via the overlaps between the wavefunctions (j4.6p with p arbitrary and those with 
Pc = |, that eq. (I2.22P associates to a climbing scalar right after the initial singularity, leads to 
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(C.2) 



It is straightforward to verify that for short wavelength modes, given the standard asymptotic 



See [48] for a similar observation related to an inflaton undergoing a random walk. 
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forms [55] 




(C.3) 



\ak\ — )• 1, |/3fc| — )• as A; — )• oo. As expected, high-frequency modes are not sensitive to the global 
expansion, and therefore the corresponding vacua should coincide. 

Modulating the power spectrum with the factor (jC.ip for the examples plotted in fig. O 
initial states natural for a radiation-dominated epoch (p = ^) rather than for the climbing 
phase would result in the spectra of fig. [TTJ Therefore, the long-wavelength tail of the power 
spectrum would be enhanced to the point of almost nullifying the suppression induced by the 
early climbing phase, while significant long-wavelength oscillations are still present. As one 
dials up p in eq. (1021) . the effects become even larger and persist to shorter scales, although this 
would not correspond to any reasonable choice for the initial state of the Universe. All in all, if 
the Universe began in a state that is not natural for the climbing trajectory right after the initial 
singularity, the suppression of power at long wavelengths would be nullified, or even enhanced 
at the very longest wavelengths. However, such an initial state is not natural for our system, 
and we further remark that one would have arrived at similar long-wavelength modifications for 
the power spectra derived around other suitably inflating solutions (and not just the climbing 
solution), as is evident from the deviation from unity of eq. (jC.ll) for small k. This is none other 
than the usual sensitivity of CMB observables to the pre-inflationary state if inflation lasted not 
much longer than it needed tc@. 



^^In concrete terms, one can show that the factor (|C.1|) can readily introduce modulations of the power spectrum 
of order one or larger over a span of wavelengths that would have exited the horizon up to over 4-5 e-folds worth 
of inflation. We would be sensitive to this modulation if inflation had not lasted more than 63 e-folds. 
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